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Abstract 

In this paper a Time Division Multiple Access (TDMA) based Random Access (RA) channel with Successive Interference 
Cancellation (SIC) is considered for a finite user population and reliable retransmission mechanism on the basis of Contention 
Resolution Diversity Slotted ALOHA (CRDSA). A general mathematical model based on Markov Chains is derived which makes 
it possible to predict the stability regions of SIC-RA channels, the expected delays in equilibrium and the selection of parameters 
for a stable channel configuration. Furthermore the model enables the estimation of the average time before reaching instability. 
The presented model is verified against simulations and numerical results are provided for comparison of the stability of CRDSA 
versus the stability of traditional Slotted ALOHA (SA). The presented results show that CRDSA has not only a high gain over 
SA in terms of throughput but also in its stability. 

I. Introduction 

While the application of RA techniques for data transmissions is appealing in many application scenarios such as sensor 
networks, signalling or unpredictable and bursty low duty cycle user traffic, often concerns are expressed about the limitations 
in terms of spectral efficiency and the risk of RA channel instability, leading to a zero throughput and a correspondingly infinite 
transmission delay. While recently significant improvements of the spectral efficiency have been achieved by introducing SIC 
and coding techniques (e.g. [1], [2] and [3]), the stability behaviour of these new schemes has not been fully analyzed yet. In 
[4] a first analysis of the stability of CRDSA as SIC representative was done and a mathematical model for the prediction of 
the channel stability was derived, which is used as baseline for the work presented in this paper and therefore recalled in the 
following. 

The source of channel instability in a RA channel is the natural occurrence of collisions among the packet transmission and 
the presence of mechanisms which attempt the retransmission of lost packets. In principle collided packets could be simply 
discarded, but doing so would adversely affect the Quality of Service (QoS) experienced by the user or may be entirely 
unacceptable for critical signalling information, such as log-on messages. RA schemes thus usually attempt to retransmit the 
lost packets, either until they are successfully received or until a maximum number of retransmissions has been reached. In 
order to make retransmissions possible, the users need to receive feedback whether their transmission attempt was successful, 
e.g., by means of acknowledgements. The instantaneous throughput of the RA channel S{G) is then dependent on the total 
load G, being the sum of the load due to new transmissions Gp and the load due to retransmissions Gb- In this sense the 
RA channel forms a feedback loop as is illustrated in Fig. 1. It is an inherent property of closed-loop feedback systems, that 
the feedback can lead to amplifying self-excitation. Here this results in an increase of the overall load due to the additional 
retransmissions. 

The throughput curves of ALOHA, SA [5], and CRDSA [1] all have in common that for increasing load G the throughput 
S{G) first increases until reaching a maximum throughput Smax- For further increasing load, S{G) decreases again and 
asymptotically approaches zero. 

If due to retransmission attempts of lost packets the total load exceeds a critical threshold, then even more packets experience 
a collision and get lost, resulting in an even higher retransmission load. In the end the channel is driven into total saturation in 
the area of having very high load and very low throughput. To reduce this amplification effect a retransmission strategy is used, 
which shall Umit the load due to retransmissions and reduce the risk of getting more collisions (see Fig. 1). Many different 
retransmission strategies that try to achieve this goal are known from literature. In [6] the selection of the time of retransmission 
with uniform probability within a parameterizable interval t G [0, ..., K] is proposed. In [5] a strategy is described where the 
decision for a retransmission attempt is taken with a probability pr in every slot (for SA) resulting in a geometric distribution. In 
[7] the selection of the retransmission time from an interval, which grows exponentially with every collision (Binary Exponential 
Backoff), is proposed. Finally the so called splitting algorithms (see e.g., [5], [8]) iteratively split the set of collided users 
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into two sets and stabilize the system this way. Furthermore two different types of user population are distinguished, finite and 
infinite user populations. For a finite user population, every user that experienced a collision is backlogged, which means that 
he is not generating any new traffic until the collided packet has been successfully transmitted. The infinite user population on 
the other hand refers to either an infinite number of users or a finite number of users that generate new traffic independently 
of whether another retransmission is still pending or not.' While some retransmission strategies assume a visibility of the 
channel activities by all users, here we assume that every user has no instant visibility of other users activity (as is the case 
in satelUte systems with directive links and long propagation delays) and only receives feedback about the success of his 
own transmission attempt from the receiving end system. The retransmission mechanisms using a uniform and geometric 
retransmission probability have in common that the probability of retransmission is a fixed parameter and does not change 
dynamically. For the binary exponential backoff and tree splitting algorithm, the actual retransmission probability may change 
over time dependent on the situation. In the remainder of this paper the focus is on a geometrically distributed retransmission 
mechanism, since it was shown in [9] that the channel performance of SA is mainly dependent on the average retransmission 
delay and largely independent of the retransmission probability distribution. 

II. Review of SIC-based Random Access Techniques 

Over the last years the recently regained popularity of RA schemes resulted in the definition of new RA protocols. In 
particular a recent enhancement of the SA protocol, named CRDSA [1], [10], using SIC techniques over a set of slots (denoted 
frame) to improve the throughput and Packet Loss Rate (PLR) behaviour of SA, has been studied showing an impressive gain 
over SA increasing the maximum throughput from Smax,SA ~ 0.36-^ to Smax,CRDSA = 0-55;^^. Up to now however 
the consequences for the system stability of this new access scheme have not been analyzed yet. The fundamental concept 
of CRDSA is to generate a replica burst for every transmission burst within a set of Ns slots, called frame, see also Fig. 
2. While the generation of a redundant copy of a burst is similar to previous proposals such as Diversity Slotted ALOHA 
(DSA) [11], the fundamental difference here is that every burst contains a pointer to the location of its replica. In case a clean 
replica arrives, meaning that the burst could be decoded and received successfully, the channel is estimated from it and the 
interference that this burst introduces to other users is removed for all replica-burst locations. 

In the example in Fig. 2 the first burst of user 1 is received successfully since not interfered. As consequence of the SIC 
process, the interference that the replica of user 1 introduces to the second burst of user 2 is removed so that this burst of 
user 2 can be decoded in the next round. This process is then iteratively repeated. In the example in Fig. 2 all repUcas can be 
recovered this way. 

A. Characterization of the Packet Loss Rate in CRDSA 

For classical SA, the necessary condition to have a successful reception is that only a single transmission must occur in a 
timeslot, otherwise the burst is lost. Let us denote by M the total user population of the system and po the probability that 
a user attempts a transmission in a time slot, then the probability that a user successfully receives a packet gets Psucc,SA = 
Po • (1 ~ Po)^^ ^- Increasing the overall number of users A/ — > oo, the totally transmitted packets can be modeled as 
Poisson process with arrival rate A [12]. The probability for a successful transmission then results in the well known equation 
Psucc.SA = ^ 'Tp ■ e^'^ ^p, whereas Tp denotes the slot duration. This simple closed form expression is conveniently suited to 
describe the throughput surfaces, which are used for the stability investigation done e.g., by Kleinrock [13]. The preconditions 
in CRDSA are however different due to the iterative SIC process. As was shown by Liva in [2] and [14], the SIC process can 
be interpreted as an erasure decoding process in a bipartite graph, such as for Low Density Parity Check Codes (LDPC) codes 
[15]. For this purpose, every slot in a frame is represented as a sum node and every transmitted burst by a burst node. The 
edges in the graph then connect the burst nodes to the sum nodes. In [2] an expression for the average erasure probabilities 
for every iteration are derived for the asymptotic case of infinitely long frames, resulting in an upper bound of the achievable 
throughput. An expression for neither the exact nor the average erasure probabilities in a non-asymptotic case with finite 
frame lengths however can be expressed accurately by these bounds or another closed form expression. For this reason the 
stability analysis in this work relies on simulated CRDSA packet success probabilities and throughput for the case of having 
one additional replica (degree d = 2), a frame consisting of Ns = 100 slots and a limitation of the number of SIC iterations 
to Imax = 10. The presented framework is however flexible to be used as well for other configurations of CRDSA, always 
requiring only that the average throughput curve is known. 

B. Stability Definition 

The issue of stability in RA systems was already identified in the very early days of the ALOHA proposal. Abramson [16] 
and Roberts [17] both addressed this issue for plain ALOHA. After the evolution of ALOHA towards SA, many publications 
have dealt with the investigation of the stability behaviour of SA, for instance [5], [18], [12] and [13]. Stability is commonly 

'Generating a new transmission in addition to a retransmission can be also seen as two users, one retransmitting, one transmitting new data. Since the 
generation of new transmissions is not bounded, the user population can also grow to infinity. 



defined as the ability of a system to maintain equilibrium or return to the initial state after experiencing a distortion. In the 
context of RA, the term stability is used in different ways in literature. In the definition given by Abramson in [16], the 
ALOHA channel was defined instable if the average number of retransmissions becomes unbounded. Within [5] a channel was 
defined stable if the expected delay per packet is finite. Kleinrock defined in [13] a channel as stable if the SA equilibrium 
contour (i.e., throughput is equal to the channel input rate) is nontangentially intersected by the load line in exactly one place. 
In the strict mathematical definition of stability of autonomous systems, this corresponds to a sufficient condition for a global 
equilibrium point. In the terminology used by Kleinrock, a SA channel is instable if the load Une intersects the equilibrium 
contour in more than one point. In the mathematical sense also then the system can have a locally stable equilibrium point, 
so the definition of stability by Kleinrock refers to the criterion of having a single globally stable equilibrium point. In the 
remainder of this work, the definitions given in [13] are followed also here, meaning that a channel is denoted as stable if it 
has a single globally stable equilibrium point and instable otherwise. 

III. Stability in CRDSA 

Within this section, the derivation of a Markov model for a finite user population is described and the mathematical 
formulations for throughput and drift are derived, which form the core of the stability framework presented afterwards . This 
section concludes with a stability analysis for a representative CRDSA configuration. 

Let the RA channel under consideration be populated by a total of M users (finite user population). Every user resides 
either in a so called fresh (F) state or backlogged (B) state. In the beginning all M users are in state F. Every user in state 
F attempts a new transmission in the current frame with probability pq. It is further assumed that all users receive feedback 
about the success of their transmission at the end of a frame. In case the transmission attempt was successful, the user remains 
in state F. In case a packet is lost, the user enters state B. A user in state B attempts a retransmission of the lost packet with 
probability pr in the current frame. In case the retransmission is successful the user then returns to state F, otherwise the user 
remains in state B. Let X^ denote the number of users in state a £ {F, B} in frame /, then the discrete-time Markov chain can 
be fully described by either Xg or Xp, since both are connected by Xp = M — X'g. In the following X^ is chosen as the 
Markov state variable. Given the initial state Xg = and the state transition probability P{x'\x), which is the probability to 
move within one frame from backlog state x to state x', the Markov chain is then fully described. One major difference to the 
SA analysis done by Kleinrock is that the backlog state for SA can at maximum decrease by 1 user per slot (otherwise there 
would be a collision), while the backlog state Xb for CRDSA can decrease by [1, . . . , Xb] in a frame. Since no closed form 
expression for the success probability of a user in CRDSA is known in literature, the probability qd,NsJmai'''j''-') i^ introduced, 
which is the probability that out of r users who attempt a transmission in the frame exactly v users are successful. The success 
probability q is dependent on the CRDSA configuration, consisting of the repetition degree d, the number of slots in the frame 
Ns and the maximum number of iterations Imax- Here, this probability was derived numerically by simulations and averaging 
over the results for every offered load G. For sake of simplicity, the subscripts will be omitted in the following, using g(r, v). 
When changing state, let vi be a random variable denoting the number of successful transmissions in frame /, ipi a random 
variable denoting the number of fresh transmission attempts in the frame and pi the number of retransmission attempts in a 
frame. 

Let FSi denote the number of fresh users, which transmit successfully in frame I. Let FUi be the number of fresh users 
who attempted a transmission but were unsuccessful. In the same way, BSi denominates all backlogged users who attempt 
a retransmission and were successful and BUi those backlogged users, whose retransmission attempt was unsuccessful. The 
following equations (l)-(3) can be derived: ^ 

ipi = FSi + FUi, (1) 

Pi=BSi+BUu (2) 

vi=FSi+BSi. (3) 

The joint probability mass function, conditioned on state X^ = xb is then given by Eq. (4). 
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-It should be noted that idle users have no relevance here since they neither change the size of the sets Xp and Xb nor do they generate load which 
impacts the transmission perfomiance. 



With Eqs. (l)-(3) the change in number of backlogged users Aa;^ = xb,i+i — xb,i = x'g — xb = FU — BS can be easily 
reformulated into: 

V = if + Xb - x'^ (5) 

The state transition probability P{x'g\xB) can then be formulated by combining (4) and (5) into (6); 

Pix'B\xB) = 

= ^P{f,P,f + XB -x'g\xB) = 



M-xi 



E("7^)!'S(1-») 

■q{ip + p,ip + XB -x'g). (6) 

With (6) in principle the entire Markov chain can be described with all its transition probabilities. In practice the computational 
cost of computing all transition probabilities is however enormous, mainly due to the nested summations over a large range 
of possible values for ip and p. To avoid this computational complexity, the stability analysis in the following makes use of 
a drift analysis, in reminiscence of [13] and [19]. The change of the backlog state forms a differential equation, whereas the 
drift corresponds to the change of the state variable ds = -|^. For the drift analysis the change in backlog xb over time is 
analyzed in the following and the stability of the equilibrium points is computed by using the tools known from differential 
calculus. In the style of [19] and [20], the drift is here defined as the expectation of the change of the backlog state Xg frame 
by frame as given by: 

dixB) ^dB^E {X'+^ - X'b I X'b] = 

= Y,{x'b-xb)-P{x'b\xb)- (7) 

•^ B 

whereas E{.} denotes the expectation value E^\y{f{x)} = J2fi^) ' Pi^lv)- 
With (5) and (6) this can be reformulated into: 

d{xB) = ^if-v)- P{x'b I xb) = 

V-P.v 

= X! ("^ ~ ■") ■ ^('^' P' ■" I 2;b) = 

= S{$}-^{T}. (8) 

whereas $ denotes the random variable taking the values (p and T the random variable taking the values v. 
From (4) it is clear that (p is binomial distributed so: 

E{<P}^iM-XB)-po- (9) 

The second expectation value E{T} is related to the throughput S{xb) of the system by (10), i.e., the expected number of 
successful packets per slot in frame /: 

S{xb) = -^E{T} = 
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The expected throughput S{xb) can also be expressed via the average success probability Ps{x), i.e. the probability of a 
successful transmission in a frame when attempting x transmissions: 

S{xb) = [{M ~ Xb)PO + XBPr] ■ Psi{M - Xb)PO + XBPr)- (H) 



With (9), (10) and (11) the drift ds becomes: 

ds = (M - xb) -Po-Ns- S{xb). (12) 

With (12) it is now possible to fully describe the stability of the CRDSA system for the case of having a user population M, a 
probability po of fresh users generating new packets and a retransmission trial probability of pr- Intuitively, the drift represents 
the tendency of the system to change over time and gives the direction of change of the backlog size. This means that for 
positive drifts the size of the backlog tends to increase by ds (i.e., more users experience lost packets and get backlogged). 
For negative drifts, the length of backlog decreases, which means that backlogged users successfully retransmit and get fresh 
again. A drift of corresponds to an equilibrium point, which may be locally stable or instable. Fig. 3 shows the dynamics 
of the channel with the drift-backlog surface for the scenario M = 500, Pr = 0.78 and for varying pq. The surface can be 
classified into three different areas: In the first area for < po < 0.01 the drift-backlog-surface does not intersect the zero-drift 
plane and gets the tangent plane for pq = 0.01. In the second area for 0.01 < po < 0.11 the drift-backlog-surface intersects 
the zero-drift plane in three equilibrium points. In the third area for 0.11 < po < 1 the drift-backlog-surface intersects the 
zero-drift plane in a single equilibrium point which is located at the saturation point where all or almost all Af users are 
backlogged. 

Fig. 4 shows the backlog drift of the three areas for three representative values of po, i.e. the intersection of the drift-backlog 
surface from Fig. 3 with the planes po = {0.01, 0.04, 0.12}. 

As can be seen here, the drift for the stable configuration (po < 0.01) is always negative independent of the backlog state xb 
and approaches asymptotically a drift ds = 0, which means that the system always shows the tendency to lower the current 
backlog state until reaching the initial state. There is thus only one equilibrium point (globally stable) close to the initial state. 

For the instable configuration (0.01 < po = 0.04 < 0.11) it can be seen that after the initial equilibrium point (locally stable) 
and the following area of negative drift (up to xb ~ 209) a second, locally instable equilibrium point is reached at xb = 209. 
When reaching this point the system can either fall back into the negative drift region for xb < 209 or enter the region of 
positive drift xb > 209. In the latter case the positive drift means that any movement to a higher backlog state (which is a 
consequence of the positive drift) results in an accelerated increase in number of backlogged users. This behavior then persists 
until reaching the third and final equilibrium point (locally stable) at xb ~ 500 = M. In this third equilibrium point now all or 
almost all M users are backlogged and the system has reached the point of maximum load and minimum throughput. In case 
some of the users get unbacklogged, the drift is anyway positive and drives the channel back into the saturation. While there is 
a low probability that the channel returns in the high throughput region, the probability is fairly small and it can be expected 
that a very long time passes before this happens. For the overloaded configuration 0.11 < po = 0.12 < 1.0 the drift-backlog 
surface intersects the zero drift plane only once at the saturation point where all or almost all M users are backlogged. Since 
there are no equilibria before and the drift is always positive, it can be expected that the system moves straight towards the 
saturation point after being started. An instable system may remain for some time in the desirable high throughput region of the 
operating point before getting instable and entering the low throughput region around the saturation point. For an overloaded 
configuration the channel moves directly to the saturation point. 

From this observation the conclusion can be drawn that for a given system configuration ft — {d, Ns,Imax} the maximum 
traffic generation probability po for which the system is still always stable is the one resulting in a drift contour which intersects 
the straight line d^ = at most once (Fig. 4). The resulting single equilibrium point is then locally and globally stable. For 
all other cases the channel is instable (e.g., instable configuration with po = 0.04 in Fig. 4), meaning that earlier or later 
the backlog will increase into the total saturation point. If the single equilibrium point coincides with the saturation point, 
the system is overloaded. While it is - mathematically speaking - also stable in this scenario (locally and globally stable 
equilibrium) it is in total saturation with very low throughput and high delay. Since an operating point in this region is not 
viable for a communication system, in the remainder of this paper stable refers only to having one equilibrium point in the 
high throughput region. 

With this framework, it is now possible to predict the stability of a channel with a certain set of parameters or to derive 
a set of parameters for which the channel is guaranteed to be stable. In [4] the validity of this model was verified against 
simulations in different scenarios. 

IV. Average Delay 

From the stability model defined in the previous section it can be observed that the stability of a system with fixed pq 
benefits from a reduction of the retransmission probability pr- Or in other words a configuration which is instable can always 
be stabilized by decreasing p,.. This comes however at the cost of a higher delay since reducing pr means increasing the average 
time before attempting a retransmission. On the one hand a low average delay (i.e. requiring pr to be as high as possible) is 
important for achieving a good QoS perception for the user On the other hand remaining stable is important for user satisfaction 
as well, since an instable system will be driven into total saturation with asymptotically zero throughput and infinite delay. 
For a stable configuration however it is beneficial to have a p,. as small as possible. The retransmission probability pr is thus 
a design parameter which can be optimized to achieve a delay as low as possible while being selected high enough to ensure 



a stable system operation. For this reason it is important to derive an analytical framework that makes it possible to compute 
the expected delay for a given system configuration in order to find the optimum choice for the design parameter p.,,, e.g., to 
minimize the delay while remaining stable, but also for optimizing the maximum allowable packet generation probability po 
or the maximum allowable user population M which is treated in section V. While the stability model derived in the previous 
sections provides the mathematical framework to derive the overall set of parameters for which the RA is stable, this section 
deals with the computation of the expected delay for any set of parameters. 

The analysis of SA in [13] followed the fundamental principles of Markov theory and derives the expected delay Di, via 
Little's theorem. According to this well known theorem, the average number of packets in a queuing system in stable conditions 
is the product of the packet arrival rate and the average dwell time in the queue. Applied to the stability analysis, the expected 
dwell time in the queue corresponds to the transmission delay Di, of every packet (i.e. the time the packet remains in the 
channel until it is successfully received). The average number of packets in the channel is given by the expected backlog length 
N since every backlogged user has one pending transmission. In equilibrium the traffic arrival rate is equal to the serving rate, 
or in other words the channel throughput 5*0 is equal to the offered load Gq. Following this analogy, the expected delay in a 
random access channel computes with Little's theorem to: 

TV 
Db = ^ (13) 
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and Pn is the probability of being in state n. Similarly the expected backlog length N can be computed as: 
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As it has been shown in [13] by numerical simulations, the values for Sout and N can be closely approximated by the 
equilibrium point throughput 5*0 and backlog state no, i.e. Sout ~ Sq and N ^ uq. With this and (13) the expected delay Dt 
gets: 

A = ^. (14) 

In order to show that the approximations of Sout and N claimed by Kleinrock for SA are also valid in the case of CRDSA, 
the theoretical expected Di, and the measured delay £)*^™ have been compared for a representative CRDSA configuration 

* = {M,po,Pr,d, TVs, /,naj = {200, 0.9, g^, 2, 100, 10}. 

The channel is stable in this configuration with an equilibrium point at no = 149.06 and an average throughput of 5*0 = 0.46. 
With (14) the expected delay gets Df, = 324.04 slots = 3.2404 frames. The average delay obtained by simulations is £)*"" = 
328.64 slots = 3.2864 frames, which is fairly close to Dt and thus confirms firstly that the approximation for Sout = So and 
N = no are also valid in the case of CRDSA, and secondly that the presented framework is suitable to estimate the average 
delay for a given channel configuration ^. 

V. Stability Comparison of SA and CRDSA for Stable Channels 

With the ability to compute the expected delay for a given configuration ^, the stability of CRDSA can now be compared to 
the stability of SA. The stability of SA was deeply investigated in [13]. The comparison of the two stabilities is of particular 
interest since CRDSA offers much higher throughput rates, also for higher offered traffic loads but the question arises whether 
this gain comes at the cost of lower stability, or not. For comparing the two RA schemes it needs to be ensured that the 
conditions are comparable. For the stability and performance of the RA schemes a tradeoff exists between the total user 
population-* M, probability of traffic generation po for unbacklogged users (user activity) and retransmission probability pr for 
backlogged users. As it was explained earlier, the RA channel for a finite user population can always be stabilized by choosing 
a low enough retransmission probability pr- The selection of pr on the other hand impacts the delay, as shown before, e.g., 
a lower value of pr will have a positive impact on the stability of the system but results in longer delays. For comparing the 
SA and CRDSA stability the following optimization criteria can be chosen now: 

1) Minimize the average delay Dt for fixed user population M and fixed traffic generation probability po 

2) Maximize the size of the user population M for fixed po and average delay Db 

3) Maximize the supported traffic generation probability po for fixed M and Di, 

'Here only finite user population scenarios are considered 



In the following analysis of these three criteria, CRDSA configurations are specified by the set $ whereas the SA configurations 
are denoted by the set S = {M,pQ,pr}. It should be noted that the traffic generation probability pq refers to the probability of 
generating a packet in a transmission frame for CRDSA whereas it refers to the probability of generating a packet in a time 
slot for SA. To ensure a fair comparison among the two, i.e., having the same overall traffic generation, pq for SA is chosen 
to (15) in the following. 

„CRDSA 

P'o^ = '-^^ (15) 

A. Comparison of Achievable Delay Df, 

The design parameter pr impacts the stability of the system as well as the resulting average delay. For a given user population 
M and traffic generation probability po this forms an optimization problem of selecting the optimum p* which is low enough 
to guarantee a stable operation of the channel, while it should be as high as possible at the same time to provide a low average 
delay. This optimization problem can be formulated for CRDSA in the following way: 

p* = argmin Db{'^,pr). (16) 

Pre[o,...a] 

Fig. 5 illustrates this optimization problem. As it can be seen, the argument resulting in the minimum achievable delay for 
CRDSA and configuration *' = {M = 400, po ^ 0.263, Pr,d = 2,Ns = I00,lrr,ax = 10} gets; 

Pr,CRDSA = argmin Db{'i'',Pr) = 5 • 10"^ 

P.e[o...i] 

and the resulting minimum average delay computes to: 

D'^^]{j)r = P*. CRDSA = 5 • IQ-^) = 3.68 frames = 368 slots. 

For SA and comparable configuration S' = {M = 400, po = 2.63 ■ 10~^, Pr} the optimization for the minimum achievable 
delay results in the optimum retransmission probability p* ^^ = 2.5 ■ 10"'^ and an average delay of: 

D^siPr = Pr,sA = 2-5 • 10"^) = 707 slots = 7.07 frames 

While naturally the gain in terms of delay of the two different schemes changes with the other configuration parameters in 
^P and S, the results above show that in the given configuration CRDSA can save 48% of the delay compared to SA in same 
conditions and guaranteeing a stable channel. 

B. Comparison of Supported User Population 

The second optimization criterion is to determine the maximum user population M which can be supported with the same 
average transmission delay Df, while guaranteeing the stability of the channel. Finding the maximum user population M* for 
achieving an average delay I?" forms an implicit optimization problem of Di,{pr, M) with side condition (17): 

g{Pr, M) =. DbiPr, M) - Dl = 0. (17) 

The solution of this optimization problem can be easily found with a Lagrange auxiliary function (18): 

L{pr,M,X) = Db{pr,M) + X[Di,iPr,M) - i?0]. (18) 

The maximum supported user population M* for retransmission probabiUty p* is then given simply by solving the set of 
equations: 

dpr \ 

dL{pr,M.\) I 

ax / 

resulting in the locus of tuples {M*{D");p*{D")) shown in Fig. 6 for SA with S = {M, 2.63 • 10"^, p^} and for different 
values of D'^. 

As it can be seen, the maximum user population, which can be supported at a maximum delay of Z?" = 300 slots gets 
Mgj^ == 250 for an optimum p*sa = 6 • 10"^. 

Fig. 7 shows the solution of the same optimization problem for CRDSA. The traffic generation probability was set to the 
equivalent value p'-^^^SA _ 2 53 . \q-^ \^ order to get the comparable traffic generation probability as in SA, resulting in the 
configuration '^ = {M, 2.63 • 10~^,pr, 2, 100, 10}. As it can be seen here, the joint optimization results in A^^adsa ~ ^63 
for a retransmission probability p* qrdsa ~ 0.06. 

The comparison for this configuration shows that CRDSA can support 45% more users than SA while achieving the same 
average delay and being also guaranteed stable. 



C. Comparison of Supported Traffic Generation Probability 

For the third optimization criterion, the user population is fixed together with the average delay to be achieved while 
guaranteeing at the same time that the channel remains stable. The optimization problem here is very similar to the previous 
one in section V-B and consists in finding the retransmission probability p* for which the the traffic generation probability Pq 
is maximized for given user population M. Defining the Lagrange auxiliary function: (19) 

L{pr,P0, A) - DbiPr,Po) + X[Dbipr,Po) - D'^] (19) 

and solving the set of equations given by: 

VL(p^,Po,A) =0 

provides the locus of optimum tuples {pQ,p*) shown in Fig. 8 for SA and in Fig. 9 for CRDSA and different _D". 

As it can be seen by comparing SA and CRDSA for e.g. Z?° = 350 slots, the traffic generation probability supported by 
CRDSA is with p^ crdsa = 0-84 a factor 2.8 higher than the one for SA with p* g^ = 3 ■ 10^^. CRDSA thus allows users 
to generate traffic with a 2.8 times higher traffic generation probability than SA. 

VI. Average time before failure for CRDSA 

The previous sections were focused on the investigation of stable channels. In many application scenarios instability may be 
acceptable if the time before getting instable is only sufficiently high. In stable channel conditions, the performance comparison 
of RA schemes could be done by comparing the minimum achievable delay Di,, the maximum number of supported users 
M or the maximum traffic generation probability po- But for an instable channel configuration, these criteria do not apply 
anymore. In an instable channel the operating point will sooner or later reach the locally stable but undesired equilibrium point 
Us in the low throughput region, which can also coincide with the total saturation point, where all users are backlogged and 
the delay Df, grows to infinity. Also the maximum user population A/ or maximum traffic generation probability po are no 
suitable measures. In an instable channel, M and/or po can grow arbitrarily while the channel will always remain instable. 
What changes is the time to reach n,, which will be shorter with growing AI and po. With this in mind, the average time 
before the channel enters the instable region for the first time can be used as a suitable measure for comparing the behaviour 
of different RA schemes in instability. 

Once the undesired operating point rig in the high load/low throughput area is entered, the channel may remain there potentially 
for a very long time (unless it is being reset). There is only a very small, but non-zero, probability to get out of this undesired 
operating point which depends on the configuration $. As explained and shown in section III an instable system has three 
operating points, two of them locally stable and one locally instable. Among the two locally stable ones, one resides in the 
low load area (desired operating point) whereas the other one resides in the high-load/low throughput region (i.e. high number 
or all users backlogged). When in the locally instable operating point, the system has a chance to fall back into the desired 
region but the same chance to enter the undesired region, ending up in the low throughput operating point. In the SA analysis 
done by Kleinrock [13], this instable operating point is also denoted as the critical system state ric- A measure for comparing 
the stability of different RA channels is then the average time before the critical state ric is reached for the first time, assuming 
further that the system will fall into the low throughput region, once Uc is reached. In the Markov chain representation, the 
state ric + 1 = n„ is modeled as an absorbing state in order to simplify the analytical analysis. It should be noted that this is 
clearly only a model since in a real system the probability of leaving the high backlog state is non-zero, while it gets zero 
when using an absorbing state. In this work the focus is only to derive the time until the system is entering the instable state 
for the first time without looking at the time until it would leave the instability region again. The average first entry time Ti 
into state i can be expressed recursively by 

Tic 

T,^l + Y,p,,-T, (20) 

whereas pij denotes the state transition probabiUty from state i to state j 
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For the computation of the First Entry Times (PET), it is now of interest to know the average time until reaching the critical 
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state Tic for the first time when starting from the initial state x° = 0, i.e., T'q\ The recursive formulation in (20) yields a set 



of linear equations: 
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which can be expressed in matrix vector notation by Eq. (22) 

t = e + P t 
whereas e is the unity vector The vector of interest with all the FETs for every state i gets then 



(22) 



(23) 



with I being the identity matrix and e the unity vector For the stability measure of the channel the entry of interest is the first 
entry in t which represents the time Tq^ to reach Uc starting from the initial state xb = 0. 

A. Validation of the Model 

In order to illustrate the validity of the derived model, the markov state transition matrix and the FET is computed for a 
representative example here. The considered scenario is ^Pq = {300, 0.19, 0.7, 2, 100, 10} which was chosen to result in an 
instable CRDSA channel with equilibrium points at no = 20 (locally stable desired operating point), Uc ~ 40 (locally instable 
equilibrium) and Ug = 247 (locally stable undesired operating point). Since n„ = ric + 1 is an absorbing state it is sufficient 
to compute the Markov state transition probabilities in the range of states from [0 . . . n„]. Fig. 10 and 11 show the computed 
and simulated Markov state transition matrix P for $o- 

As it can be seen from the two graphs, the transition probabilities resulting from the simulations match very well with the 
ones derived by numerical computation. 

By solving the set of linear equations from (23), the computed FET time in this example results in FETcomp = 14.57 frames, 
where the average simulative FET reaches FETgim — 13.59 frames, which is very close to the expected FET derived by 
computation. 



VII. FET COMPARISON BETWEEN CRDSA AND SA 

For a fair comparison between the FET of CRDSA and SA, configurations need to be selected which have the same initial 
conditions, i.e., the same user population M and traffic generation probability po- The average delay D[, cannot be used here 
since the average delay Di, is infinite for an instable channel by definition as the channel will enter the saturation point with 
close to zero throughput. 

A difficulty here consists in the fact that for an instable CRDSA configuration, the SA channel is getting overloaded. 

On the other hand an instable SA configuration for which three equilibria exist results in a CRDSA configuration which is 
stable so no FET can be computed. 

Fig. 12 shows the FET times for different configurations of CRDSA and SA for a user population of M = 500. 

As it can be seen, the FET for CRDSA and pC^R^SA _ iq-i jg ^p [q ^ factor 20 higher than for the equivalent SA 
configuration. It should be noted that SA is already in overload for a Pq"^ = 10"'^. For this reason the SA FET curve does 
not show the time until reaching the critical state ric but the time until reaching the saturation point Us instead. 

Fig. 12 furthermore shows the FET curve for CRDSA until reaching the saturation point n,,. The FETs for this curve are 
slightly higher than for Uc as could be expected. This result also confirms the assumption to compute the FET times by 
modeling ric as an absorbing state instead of computing the full Markov chain up to Ug, since once ric is reached also Us is 
reached very fast. The SA curve in Fig. 12 could arise the impression that the FET curve is flat and has a qualitatively different 
shape than the CRDSA curve. This is actually not the case and the Tq" for p,. — 0.5 is indeed higher than for pr = 0.9 with 
a value of T^'lpr = 0.5) = 6809.04 slots and T^"" {p,- = 0.9) = 6798.45 slots. In this configuration the SA channel is already 
so overloaded that also a large decrease of the retransmission probability pr does not have a significant impact anymore. Once 
the traffic generation probability p^^ is lowered, the impact of pr gets more visible as it can be seen from the last curve in 
Fig. 12, which was computed for Pq^ = 3.57 • lO^"'. Also for this lower Pq^ the Tq" is much lower than for CRDSA with a 



higher Pq 



CRDSA 



10 ^, showing that for an instable configuration CRDSA is remaining stable much longer than SA. 



VIII. Summary and Conclusions 

In this paper, a theoretical model for the stability of CRDSA as representant for SIC RA schemes was developed. With 
this model it is possible to draw qualitative and quantitative conclusions about the stability of the communication channel. 
The presented framework enables the estimation of the average delays experienced in stable channel configurations. The stable 
CRDSA and SA RA channels were optimized for achieving a minimum delay, maximizing the user population while achieving 



10 

a delay target or deriving the maximum traffic generation probability for a given user population and delay target for which the 
channel is stable. Numerical results were presented which allow a direct comparison of the performance of CRDSA and SA. 
These results have shown that CRDSA does not only provide a higher throughput and lower PLR than SA but is also capable 
to achieve lower delays and higher user population and traffic generation probabilities than SA while being stable. Finally the 
stability framework was extended towards instable channel configurations of CRDSA and makes possible to predict the average 
time before reaching instability. The derived model for CRDSA was validated against simulations and the stability behaviour 
of CRDSA was compared to the one of SA for instable channels. Also here CRDSA showed a much better performance by 
reaching way higher average times before failure than SA. Besides the analysis of the stability behaviour of a channel, the 
presented framework enables the computation of the optimum design parameters, in particular p,. for which the channel either 
remains guaranteed stable while minimizing the average delay or the pr for an instable channel, which results in the desired 
FET time. 
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Fig. I . Feedback loop of a RA channel. 



Fig. 2. SIC principle of CRDSA. Numbers at an'ows indicate tlie step in tlie iteration process. 
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Fig. 3. Drift backlog surface for M = 500, pr = 0.78. 
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Fig. 4. Backlog drift for a stable channel with M = 500, p,- = 0.78 and po = 0.01, an instable channel with M = 500, pr = 0.78 and po = 0.04 and 
an overloaded channel with M = 500, p,- = 0.78 and po = 0.12. 
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Fig. 5. Comparison of minimum achievable delay for SA and CRDSA and "I-' = {M = 400, po = 0.263, pr,d = 2,7Vg = 100, /^aa; = 10} and 
S' = {M = 400, PO = 2.63 ■ 10-3, p^}. 




Fig. 6. Locus for Di,{M,pr) = D'j^ with Djj G [180 . . . 340] slots for SA and locus of max. supported user population A/(D[J). 
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Fig. 7. Locus for Db{M,pr) = -DP with D" G [200 . . . 340] slots for CRDSA aiid locus of max. supported user population Af (D"). 
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Fig. 8. Joint optimization of po and pr to achieve D^ slots for SA with H = {250,poiPr}. 
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Fig. 10. Mai'kov state transition probabilities pij for 'I'o. 
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Fig. 11. Simulated Markov state transition probabilities pij for ^o (1000 runs). 




Fig. 12. FET times for CRDSA and SA. 



